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Abstract 


An  integral  representation  for  the  particle  velocity  in  terms  of  a 

Green's  function  and  certain  linear  combinations  of  the  inelastic  strain 

■,  <  .  ">  '  -  » 

rates  has  been  used  in  this ^papet-,  both  for  a  numerical  method  to  compute 
full  field  solutions  and  to  develop  unequivocal  asymptotic  expressions  for 
the  near-tip  fields.  Specific  results  have  been  obtained  for  a  stationary 
crack  in  a  solid  whose  constitutive  behavior  is  represented  by  the  Bodner- 
Partom  model.  It  is  shown  that  the  leading  term  of  the  near-tip  particle 
velocity  is  of  order  r  ,  and  the  higher-order  terms  are  of  the  forms  r  logr 
and  r.  Expressions  ^have  been  derived  for  the  angular  variations  and  for  the 
multiplying  time-dependent  intensity  factors.  The  r  logr  term  is  absent 
for  the  Mode- II  case.  Two  questions ,, have-been > addressed  in  further  detail: 
the  dependence  of  the  multiplying  terms  on  time  and  the  importance  of  the 
higher-order  terms.  The  numerical  results  show  a  stress  intensity  factor 

which  decays  with  time.  At  a  small  distance  from  the  crack  tip  the  numerically 

.  u  •>  i 

computed  normalized  opening  stress  ^22^xi,®,t^T22^xl'^,^^  ^a8  been /comPare<* 
with  a  one-term  asymptotic  representation,  i.e.,  with  K^(t)/KI(0) . '>The  two 
curves  diverge  at  very  small  values  of  time.  "The/inclusion  of  a  second  term 
in  the  asymptotic  expression  for  the  stress  gives  very  acceptable  agreement 


as  time  increases. 


1.  Introduction 


A  clear  understanding  of  the  mathematical  structure  of  the  fields  of  stress 
and  deformation  near  a  stationary  crack  tip  is  an  essential  prerequisite  to 
an  investigation  of  conditions  for  the  onset  of  crack  propagation.  In  a 
visco-plastic  solid  the  fields  in  the  immediate  vicinity  of  a  crack  tip  can 
be  very  intricate.  In  this  paper  near-tip  visco-plastic  fields  have  been 
investigated  in  detail  for  a  class  of  constitutive  models  for  which  the 
inelastic  strain  rate  approaches  a  bounded  value  as  the  stresses  increase 
beyond  bounds.  The  Bodner-Partom  model  [1],  which  has  been  shown  to  describe 
the  mechanical  response  of  a  wide  class  of  materials,  displays  this  particular 
behavior  of  the  inelastic  strain  rate. 

Fields  near  a  crack  tip  are  usually  analyzed  asymptotically  by  assuming 
a  near-tip  field  of  the  general  form  rPf(0).  When  this  expression  is  sub¬ 
stituted  into  the  governing  equations  and  into  appropriate  boundary  con¬ 
ditions,  and  terms  of  equal  order  in  r  are  collected,  there  results  a 
nonlinear  eigenvalue  problem,  which  yields  both  p  and  f(9).  For  visco¬ 
elastic  solids  this  method  was  used  by  Riedel  [2]  and  Riedel  and  Rice  [3], 
who  investigated  the  near  tip  fields  of  stationary  cracks  for  Mode  III  and 
Mode  I,  respectively.  They  considered  a  viscoelastic  constitutive  equation 
of  the  power  law  type.  For  the  Mode  I  case  their  asymptotic  solution 
displays  the  well-known  HRR  singularity  if  the  power  is  greater  than  one. 

When  the  power  is  less  than  one,  the  solution  becomes  square-root  singular. 
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Hui  &  Riedel  [4]  considered  quasistatically  moving  cracks  for  both  Mode  I 
and  Mode  III  conditions.  When  the  power  is  greater  than  3,  their  singularity 
is  different  from  both  the  HRR  and  the  square-root  singular  forms,  while  for 
a  power  less  than  3  it  reduces  to  a  square-root  singularity.  Lo  [5] 
followed  the  Hui-Riedel  approach,  but  he  took  into  account  the  effect  of 
inertia. 

The  implications  of  square-root  singular  as  crack-tip  fields  for  visco¬ 
plastic  solids  have  been  explored  by  Hart  [6]  and  Freund  and  Hutchinson  [7]. 

For  a  Mode  III  quasi-static  steady-state  analysis.  Hart  [6]  related  the  stress 
intensity  factor  at  the  tip  to  the  corresponding  far-field  quantity.  He  also 
discussed  transient  cases  on  the  basis  of  certain  assumptions.  Freund  and 
Hutchinson  [7]  developed  a  similar  relation  between  the  tip  and  the  far-field 
stress  intensity  factors  for  the  dynamic  steady-state  Mode  I  case,  using  an 
essentially  linear  relation  between  the  stress  and  the  inelastic  strain  rate. 

In  this  paper  a  more  fundamental  approach  towards  the  computation  of  near¬ 
tip  fields  is  presented,  in  the  sense  that  asymptotic  forms  are  not  assumed 
a-priori.  Instead,  the  particle  velocity  field  is  expressed  by  a  representation 
integral  over  the  complete  domain  of  the  cracked  body.  This  integral  involves 
an  elastic  Green’s  function  and  linear  combinations  of  the  inelastic  strain 
rates.  Since  the  inelastic  strain  rates  are  bounded  at  the  crack  tip,  the 
integral  can  be  expanded  in  terms  of  the  distance  to  the  crack  tip,  to  yield 
an  unequivocal  expression  for  the  near-tip  field.  Details  have  been  worked 

out  for  a  stationary  crack  tip.  It  is  shown  that  the  leading  term  of  the 

u 

particle  velocity  is  of  order  r  ,  and  the  next  terms  are  of  the  forms  r  logr 
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and  r.  Expressions  have  been  derived  for  the  angular  variations  and  for 
the  multiplying  time-dependent  intensity  factors.  The  r  logr  term  is  absent 
for  the  Mode-Ill  case.  Two  points  are  of  Interest:  the  dependence  of  the 
multiplying  terms  on  time  and  the  importance  of  the  higher-order  terms. 

The  representation  of  the  particle  velocity  in  terms  of  an  appropriate 
Green's  function  and  the  time-rates-of-change  of  the  inelastic  strains  has 
not  only  been  used  to  obtain  near-tip  asymptotic  expansions  of  relevant  field 
variables,  but  also  to  develop  numerical  solutions  by  the  use  of  the  boundary 
integral  equation  method,  in  conjunction  with  an  iteration  technique.  The 
numerical  calculations  have  been  carried  out  under  the  assumption  of  small- 
scale  yielding.  Thus,  the  equations  of  the  Bodner-Partom  model  have  been 
used  in  a  region  around  the  crack  tip.  The  conditions  on  the  boundary  of 
this  region  are  provided  by  the  singular  term  of  a  corresponding  elastic 
solution.  The  results  show  a  stress  intensity  factor  which  decays  with  time. 
At  a  small  distance  from  the  tip  the  numerically  computed  normalized  opening 
stress  T22^xi,®,t^T22^xl’^’^  has  been  compared  with  a  one-term  asymptotic 
representafation,  i.e.,  with  K^(t)/K^(0).  The  two  curves  diverge  at  very 
small  values  of  time.  The  inclusion  of  a  second  term  in  the  asymptotic 
expression  for  the  stress  gives  very  acceptable  agreement  as  time  increases. 

As  early  as  1971  the  boundary  integral  equation  method  (BIEM) ,  also 
called  the  boundary  element  method  (BEM) ,  was  applied  by  Swedlow  and  Cruse 
[8]  for  elasto-plastic  material  behavior.  In  recent  years  the  method  has 
been  used  extensively  for  the  numerical  analyses  of  viscoelastic/plastic 
solids  of  various  shapes  including  bodies  containing  cracks.  Recent  work 
was  summarized  by  Mukherjee  [9] . 
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2.  Formulation 

The  states  of  deformation  that  will  be  considered  are  two-dimensional, 

either  plane  or  antiplane  strain.  The  body  contains  a  straight  crack.  At 

one  of  the  crack  tips,  say  A^,  the  applied  loads  give  rise  to  elastic  stress 

£  £  £ 

intensity  factors  K^Ct),  K^(t)  and  K^^Ct),  ^  the  material  were  linearly 
elastic.  We  will,  however,  consider  a  material  that  displays  inelastic 
behavior  as  the  stress  level  increases.  Hence  there  will  be  a  domain  of 
inelastic  deformation  near  the  crack  tip.  It  is,  however,  assumed  that  this 
domain  is  small  as  compared  to  the  crack  length,  and  that  linear  elasticity 
applies  outside  the  near-tip  domain  of  inelasticity. 

To  analyze  the  field  right  near  the  tip  A^,  we  follow  the  arguments  that 
have  been  introduced  for  small-scale  yielding.  Thus,  we  magnify  the  geometrical 
scale  near  A^,  so  that  the  geometry  becomes  one  of  a  semi-infinite  crack  in  a 
full  space,  as  shown  in  Fig.  1.  A  cartesian  coordinate  system  is  centered 
at  the  crack  tip.  For  the  plane  strain  case,  possible  boundary  conditions 
on  the  bounding  surfaces  of  the  body  are  replaced  by  the  asymptotic  condition 

[KI(t)?I(0)  +  KII(t)~II(0)]  (2>; 

as  r  •*  <*>,  where  p  is  the  shear  modulus,  and  (r,9)  are  polar  coordinates 
centered  at  the  tip  in  the  magnified  configuration.  Note  that  the  definition 
of  the  stress  intensity  factors  differs  from  the  conventional  one  by  the 
factor  (2tt)  .  The  functions  f^  and  f^  give  the  angular  variations  of  the 
cartesian  components  of  the  elastic  near  tip  displacements  for  modes  I  and  II, 


fj.ce) 

<n«> 


/  <K-I> 
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i  i 

ain0/2  -  j  sin30/2  ' 

sin0/2  +  j  sin30/2  \ 
cos0/2  -  coa30/2  ) 


> 


where 

A+3p 

w-  at  ~ 

A+U 

is  a  constant  in  terms  o£  the  Lamp's  constants  A  and  u  . 

We  now  proceed  to  the  governing  equation  for  the  problem.  The  strain 


rate 

3e  t  _ 

•  6 

is  usually  decomposed  into  the  elastic  part  g  and  the  inelastic  part 
.i 

e  as 


The  elastic  part  is  related  to  the  stress  rate  x  by  Hooke's  law, 
x  ■  Altree  +  2pge  , 


where  1  is  the  unit  tensor  and  "tr"  is  the  3-D  trace.  Since  the  effect  of 


inertia  is  not  considered,  the  equilibrium  equation  applies: 


V*  T  «*  0  . 


An  equivalent  form  of  (2.7)  in  terms  of  the  displacement  rate  u  follows 

from  (2.4),  (2.5),  (2.6)  and  (2.7)  as 

A*u  **  7‘a1  , 

x-  ~  - 


where 


A*U  -  +  (X+u)  77  *u  , 

V—  **  " 


and  for  convenience  we  have  introduced  the  quantity  a  by 


a1  ■  litre1  +  2yei 


Equation  (2.8)  is  the  governing  equation  for  our  problem.  This  equation, 
together  with  the  boundary  conditions 


£:*  x*n  ■  0  on  *  0“  ,  <  0  ,  (2 

E 

and  an  asymptotic  condition  of  the  form  (2.1)  with  u  and  ^  replaced  by 
•E  "  E 

u  and  Kj  II>  defines  u.  When  Kj  ^  ■  0,  the  behavior  of  the  elastic  solution 
away  from  the  crack  tip  becomes 

u  -  0(r“is).  (2 

To  complete  the  formulation,  we  supplement  these  equations  with  the  regularity 
requirements  at  the  crack  tip: 

u  "  bounded  as  r  +  0  (2. 

x  ”  o(l/r)  as  r  H  (2. 

In  the  present  problem,  we  will  also  discuss  the  mode-III  case,  where 

x  „  -  x__  -  u  “0  (u,8  -  1,2).  We  prefer  to  use  the  stress  function  <j> 
a0  33  ct 

because  it  makes  the  stress  computation  easier.  The  stress  function  defines 
the  stress  components 


Finally,  the  asymptotic  condition  is 

$  -  -2r%  K^n(t)  cos|e  as  r 
When  K^.^(t)  *  0»  (2.21)  is  replaced  by 
•  — i< 

<)>  -  0(r  )  as  r  +  ®. 


In  this  case,  the  constant  in  (2.19)  is  zero  because  $(0)  ■  0,  as  follows 
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.  Representation  .  of  the  Velocity  Field 


The  term  on  the  rigte-hand  side  of  the  displacement-rate  equation  (2.8) 
is  of  the  nature  of  a  body-force  distribution.  This  implies  that  the  solution 
to  Eq.(2.8)  can  be  expressed  by  the  use  of  a  Green's  function.  Representations 
in  terms  of  Green’s  functions  are  at  the  basis  of  the  boundary  integral  equation 
method  (BIEM) ,  and  they  have  therefore  become  quite  well  known.  For  rate-type 
inelastic  analysis  the  appropriate  representations  were  first  developed  by 
Swedlow  and  Cruse  [  8 ] .  This  Section  gives  a  brief  derivation  of  such  repre¬ 
sentations  for  the  purpose  of  making  the  paper  self-contained.  For  further 
detail  we  refer  to  Mukherjee  [  9 ]  and  the  papers  cited  there. 

For  an  unbounded  body  containing  a  semi -infinite  crack,  the  Green's 
function  G(x,y)  is  defined  by 

A*S(s. y)  *  -  (3.1 

where  x  *  (x.,x-),  y  ■  (y^y.,)*  and  the  operator  A*  is  defined  by  (2.9). 

+ 

On  the  faces  of  the  crack  (x^<  0,  x^  ■  0“)  we  have 

T^Oc.y)  “  0  »  0.2 

where  T^G  stands  for  the  elastic  tractions  computed  from  G.  We  also  have 

G(x,y)  -  0(log |x| )  as  |x|  -*■  »  ,  (3.3 

G(x,y)  -  0(1)  as  jx|  -*•  0  ,  (3.4 

VxG(x,y)~  0 (1/ </| x | )  as  jx|  0  .  (3.5 

A  standard  analysis  based  on  Muskhelishvili ' s  method  yields  G(x,y)  as 


Gll  +  i  G2L 
G12  +  i  G22 


2ttu(1+k) 


1\  (z  -z  )  (z-z) 

I  - 2 - 2 -  -  Klog(/z-/z  )  (v  z  -  yz  ) 

,i/  iVt/z  (/z  Wz)2  0 

oo  , 


-<2log(/z+/z  )  -  log (/ z  + 

o  u 


'  \  yz  (/z+</z  )  /z(/z  -Wz)  /  2(z-z  ) 

o  o  o  o  u 

where  <  is  defined  by  (2.3),  and 

z  -  +  ix2  ,  zq  =  y  +  iy  .  (3 

Equation  (3.6)  agrees  with  the  expressions  given  by  Erdogan  [10], 

A  suitable  representation  of  the  solution  to  (2.8)  is  obtained  by  the  use 
of  Green's  identity: 

I  (A  v)  *w  dA  -I  (A*w)*v  dA  =  I  (T  v) -w  dc  -  /  (T  w) *v  dc  , 

\  e  x~  '  x  \  e  x~  '  x  3A  x~  ~  aA  x~  ' 

which  holds  for  an  arbitrary  pair  (v,w),  and  where  A  is  the  domain  bounded 

~  R )  £ 

by  circles  of  radii  e  and  R  (0  <  e  <  R)  centered  at  the  tip,  and  by  the  crack 
faces,  as  shown  in  Fig.  2. 

We  now  substitute  u(x)  for  v(x)  and  G(x,y)  for  w(x)  in  (3.8).  Subsequent 
use  of  the  divergence  theorem  and  the  regularity  conditions  (2.13),  (2.14), 
(3.4)  and  (3.5),  followed  by  letting  R  -*■  »,  e  -*■  0,  and  using  (2.12)  and  (3.3), 
gives 

u(x,t)  =  I"  V  G(y,x)  *o^(y  >t)dA^  ,  (3 

R2 

or-  in  component  form 

V?,t>  "  /  VV7,?)d*k(y’t)dAy  •  (3 

R2  k 

In  deriving  these  equations  we  have  assumed  that  the  integrals  in  (3.r' 
or  (3.10)  converge.  For  this  to  be  true,  it  is  sufficient  to  assume  that 
o’1  satisfies 


i  1 

I?  1  -  0(r  ) 


as  r  4-  0. 


(3.11) 


2  |  -  0(r  ) 


r  00 


(3.12) 


Oj_  >  -1  * 


OL^<  -1  . 


(3.13a,b) 


Although  (3.13a)  appears  unnecessarily  restrictive,  it  turns  out  that 
these  inequalities  guarantee  the  satisfaction  of  the  conditions  (2.13)  and 
(2.14)  automatically.  It  is  noted,  however,  that  (3.13)  is  just  a  sufficient 
condition . 

For  the  Mode  III  case  a  similar  procedure  gives 


(?’*)  "  (  3y  G(?»?)ejk  oJ(X.t)dAy 

p2  j 


(3.14) 


where 


g<S.Y)  “  "  77  lo8 


(/z  -  /z  )(/z  -  /Z  ) 


(/z  +  /?  )  (/?  +  /z  ) 
o  o 


(3.15) 


Hence  we  have 


a_(x,t)  “  e  3xn 


f  3  G(y,x)e  a^(y ,t)dA  . 


(3.16) 
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4.  Asymptotic  Expansions 

In  the  following  sections,  we  employ  the  potential  representation  for 

u  given  in  Section  3.  The  results  of  this  section  are  analytical,  while 

numerical  results  are  presented  in  Section  7. 

For  a  certain  class  of  constitutive  models  the  quantity  <5*,  given  by 

Eq.(2.10),  is  finite.  For  example,  a  viscoplastic  model  proposed  by  Bodner 

and  Partom  [  l]  shows  this  property.  In  this  section  we  will  obtain  the 

asymptotic  expression  for  u  for  the  case  that  lim  gi(r,6,.tr)  exists  for.  all  8. 

r+0 


In  a  conventional  asymptotic  analysis  one  seeks  near-tip  fields  of  the 

forms 


u  *  +  u2  +  u^+  - 

•  i  .  i  .  i  .  i 

a  *  ?!  +  22  +  ?3  +  - 


.i 


where  u  ^ (j  ■  1, - )  are  functions  which  are  less  singular  near  r 

with  increasing  subscript.  Since  near  r  ■  0 


g1  -  ^(O.S.t)  +  o(l) 


(4.i; 


(4.2! 


0 


(4.3 


is  bounded,  and 

V»g*  3  r^  »c)  * 

we  find  by  substituting  (4.1)  and  (4.2)  into  (2.8) 
A*u^  "  0 


(4.4, 


(4.5 


Equation  (4.5)  implies  that  is  of  the  same  form  as  for  the  linearly 
elastic  case.  The  appropriate  separation  of  variables  form  of  (4.6)  can  only 


be  guessed  at  this  stage.  One  might  try 
u2  *  r  f(9»t) 

This  expression  for  f(9,t)does,  however,  not  lead  to  a  solution  of  (4.6) 


(4.7) 


because  the  corresponding  homogeneous  problem  has  nontrivial  solutions.  In 


the  Mode  I  case,  for  example,  the  homogeneous  equation 
A*[r  f (0»t) ]  -  0  , 

with  boundary  conditions 

6  »  it!  T[r  f(0,t)]  -  0 

0-0:  i^*T[r  f(0,t)]-O  and  i2.[r  f(0,t>]  *  0  t 


(4.8) 


(4.9) 


(4.10a,b) 


where  i^  are  unit  vectors,  has  a  solution  corresponding  to  a  uniaxial  tensile 
field  in  the  x^  direction.  Hence  the  approach  based  on  (4.7)  becomes 
ambiguous  at  this  point. 

An  unequivocal  form  of  the  term  ij2  can,  however,  be  obtained  directly 

from  the  representation  given  by  (3.10).  Using  polar  coordinates  (r,9)  and 

(r  ,0  )  for  y  and  x,  we  split  (3.9)  into  two  integrals  as 

it  r 


u(r  ,0  ,t)  -  f  f  r  7G  (r,9;r  ,9  )-ai(r,0,t)drd9 
-  o  o  J  J  <  0  0 

-IT  O 

IT  <*> 

+  [  f  r  7G  (r,9;r  ,9  ) -d^r ,0 ,t)drd9 
J  J  -  ~  >  o  o  ~ 


(4.11) 


r 

c 

f 
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> 

,1 


f-: 


f£ 


I 


Here  the  gradient  operator  is  taken  with  respect  to  the  polar  coordinates 
r  and  9.  Next  VG<  and  VG  are  expanded  into  power  series  of  r.  We  can 


that  Nq  is  independent  of  9q. 


Substitution  of  (4.12)-(4.13)  into  (4.11)  yields  after  some  rearrange¬ 


ments 


$(r  ,9  )  *  L  +  +  + 

o  o  ~o  -loo  -<  o  o  ->  o  o 


(4.14) 


where 


-IT  O 


So  *  f  I  Nq(0) •ai(r,9)drd8 

7F  <*> 

// 

O 

IT 

/  </ 


(4.15) 


W 


-IT  O 
IT 


r5  N. (9,0  ) •di(r,9)drd9 
-l  o 


(4.16) 


1<(rn*9n> 

o  o 


n+1 
r  \  2 


rQ  H  ± 


sH 


39 


write 

& 

n+1 

IP 

\y 

r  <  r  :  r  VG ^  *  2 

°  "  n=0 

0 

■> 2  We-V 

(4.12) 

V  < 

\V 

$ 

oo 

n+1 

r  >  r  :  r  VG  »  2 

o  -->  , 

<r 

■>  2  W9-V 

(4.13) 

r 

n*-l 

The  coefficients  in  (4. 12) -(4. 13) 

can  be  obtained  from  (3.6). 

It  is  found 
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[2  (^-)  Mn+i(0*eo>3  -Nq(0)-N  (0,0  )(~)  } *o  (r ,0)drd0 

°  n"°  0  (4.17) 
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The  dependence  on  time  t,  shown  in  (4.11)  is  implied,  but  not  explicitly  shown 

in  (4.l4)-(4.i8)  and  in  the  sequel.  In  the  next  step  the  orders  in  r  of  the 

o 

integrals  (4.17)  and  (4.18)  are  determined. 
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if 

V-V 


m 


-.v;  \ a'.'/X'.vV.  av.v 


Since  g  (r,9)  has  an  estimate  given  by  (4.3)  it  follows  that  (4.17) 


is  of  order  r  ,  i.e., 

o 

,0  )  -  r  A(0  )  +  o(r  ) 
<  o  o  o  o  o 


(4.19) 


where 


IT 

•  r  00 

A(8o)  =■  /  ([  Z  ^3  «  +1(9,0O)3-N  (9)  -  2  N  (0,8  )}-g1(O,0)d0  (4.20) 

-it  n-0 

A  discussion  of  the  order  of  the  integral  in  (4.18)  is  given  in  Appendix  A. 

The  result  is 


f,(9)  r  log  r  +  r  B(9  )  +  o(r  ) 
-  z  o  o  o  o  o  o 


(4.21) 


where 
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IT 

/  N2(0,0Q).di(O,0)d0 

— TT 


(4.22) 


The  function  B(0q)  is  defined  by  Eq.  (A.  6).  By  combining  (4.14),  (4.19) 
and  (4.21)  it  is  concluded  that  u(r  ,8  )  is  of  the  form 

*•  O  O 


^  +  ^1(9r.)ro  +  ^^n^n  1(JS  Tn  +  t^n^n  +  0<O 
o  o  ~o  1  OO  'Zoo  ”0  ~J  O  O  o 


(4.23) 


where 

f’_(9  )  -  A(0  )  +  B(0  )  (4.24) 

-Jo  -  o  -  o 

The  constant  fQ  in  (4.23)  corresponds  to  a  rigid-body  motion;  it  will 

be  left  out  from  further  considerations.  Some  additional  manipulations 
«  • 

reduces  f, (8  )  and  f-(9  )  to  the  forms 
-1  o  ~2  o 


W  -  [KIfI(0o)  +  KIIfII(0o)y2u 


(4.25) 
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(<+l)cos0 

(<-3)sin0 


(4.26) 


where  f^(0Q)  and  £jj(9q)  are  deiined  by  (2.2a,b),  and 


KI,II(t)  "  2tt(<+1) 


1  T  —  J«  •  i 

<+l)  J  j  T  -I, II'-  ^r»0»t)drd0  » 


(4.27) 


8iry  (k+1) 


by  r  -i,n,-1(0’e’l:)'ie  ■ 


(4.28) 


In  (4.27)-(4.28) 
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,  (4.29b) 


(1  -  k)cos20  -  2cos40  , 


-2  sin40 


-  2sin40 


(1-k)cos20  +  2cos40 


(4.30a) 
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(<-3)sin20  +  2sin40, 
2(cos20-cos40) 


2(cos29  -  cos40) 
(ic-3)sin20  -  2sin40 


(4.30b) 


It  is  shown  in  Appendix  A  that  f^(9)  can  be  written  in  the  form 


~f3<V 


/■it  ..  /(<+l)cos9  \  /  sin8 

P(0;0  ).g1(O,0)d0  +  C.  °  +  C,  0 
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o  o 


(4.31) 


The  function  P(9;9q)  is  complicated.  For  the  Mode  I  case  we  find  after 


lengthy  manipulations 
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where  C  is  a  constant ,  and 
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The  symbol  x  defines  a  product  of  a  column  and  a  row  vector.  For 
example 


(4.34a,b) 


(4.34c) 
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gf 


Equation  (4.23)  implies  that  the  dominant  singularity  of  the  strain  rate 
has  exactly  the  same  form  as  in  linear  elasticity.  The  rates  of  the  stress 
intensity  factors  ^  are  given  by  (4.27).  The  second  singular  term,  is  of 
the  order  of  logr  near  the  tip.  The  magnitude  of  the  multiplying  term  is 
determined  only  by  the  distribution  of  d*  at  the  tip.  Because  of  this  last 


property,  the  expression  (4.28)  holds  true  for  any  crack  configuration, 
whereas  (4.27)  is  true  only  in  the  sense  of  the  "small-scale  yielding". 
The  logarithmic  singularity  in  the  strain  rate  does  not  give  rise  to  a 


logarithmic  stress  singularity  except  in  for  mode  I. 

For  the  mode  III  case,  we  can  use  the  same  general  approach  to  show  that 


(-1).  f»». 
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\  COSJ 
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Note  that  the  stress  rate  and  hence  the  strain  rate  do  not  have  logarithmic 
singularities.  As  for  the  mode  1  case,  we  can  relate  the  absence  of  the 
logarithmic  term  either  to  the  non  existence  of  a  non-trivial  solution  to  a 
certain  boundary  value  problem  or  to  the  lack  of  a  certain  term  in  the  far 
field  expansion  of  the  Green's  function. 

In  the  published  literature,  the  choice  of  terms  for  asymptotic  expansions 
is  generally  made  on  a  trial  and  error  basis.  The  analysis  presented  in  this 
section  gives  a  more  rational  foundation  to  asymptotic  near-tip. 


5.  -  The  Bodner-Partom  Model 

In  Ref.[l  ],  Bodner  and  Partom  proposed  a  set  of  constitutive  equations 
for  an  elastic-viscoplastic  material.  These  equations  have  the  convenient 
property  that  no  separate  specification  of  a  yield  criterion  is  required,  nor 
is  it  necessary  to  consider  loading  and  unloading  separately.  Within  the 
context  of  the  Bodner-Partom  model  both  elastic  and  inelastic  deformations  are 
present  at  all  stages  of  loading  and  unloading,  but  the  inelastic  deformations 
are  very  small  at  low  stress  levels.  The  Bodner-Partom  equations  have  been 
used  to  model  the  mechanical  behavior  of  several  metals  over  a  wide  range  of 
temperatures  and  strain  rates.  Very  satisfactory  agreement  between  theoretical 
and  experimental  results  has  generally  been  obtained.  For  a  recent  discussion 
of  the  Bodner-Partom  model  we  refer  to  Ref. [12], 

The  constitutive  equations  follow  by  supplementing  Eqs.(2.5)  and  (2.6) 


with  an  expression  for  the  rate  of  inelastic  strain.  In  a  recent  version  of 
the  Bodner-Partom  theory  this  strain  rate  is  given  by,  see  [ 12] » 


-  “  Do  7j“  exPt"<z2/3J2^ 


where  t '  is  the  stress  deviator 


In  (5.1)  and  (5.4) »  Z  ,  Z, ,  m,n  and  D  are  material  constants,  while  W  Is 

o  1  o  p 

the  work  over  the  Inelastic  strains, 

W  ■  /  x •  dt  . 

p  J  -  - 

For  uniaxial  tension  it  is  simple  to  plot  the  inelastic  strain  rate 
versus  a/Z.  For  various  values  of  n  the  curves  are  shown  in  Fig.  3. 

The  linear  plot  of  Fig. 3  shows  that  remains  very  small  at  small  stresses 
with  a  rapid  change  over  a  small  range  of  a/Z  (the  "yield"  region) ,  and  that 
an  asymptotic  limit  is  approached  as  a  becomes  large. 

It  is  of  interest  to  note  that  the  components  of  e*,  as  defined  by 
Eq.(5.1),  are  always  bounded.  Even  if  t  is  singular  (which  is  the  case  at  a 
crack  tip) ,  the  singularity  of  t  as  r  -*•  0  will  be  cancelled  by  the  singularity 
of  in  the  term  tera  exp [-(Z2/3J2)n/2]  is  bounded  whether  J ^ 

is  singular  or  not.  Hence  Eq.(5.1)  Implies  that  the  inelastic  strain  rate  is 
bounded  at  a  crack  tip.  It  then  follows  from  Eq.(2.10)  that  the  components 
of  a^  are  also  bounded  at  a  crack  tip,  and  it  can  thus  be  concluded  that  the 
asymptotic  results  of  Section  4  apply  to  the  Bodner-Partom  model. 


.  6 .  Numerical  Procedure 


A  numerical  procedure  has  been  developed  to  compute  the  stress  intensity 
factor  and  the  opening  stress  at  a  P°  sition  on  the  crack  line.  As  pointed 
out  by  Mukherjee  and  co-workers  [13],  a  numerical  evaluation  based  on  the 
representation  (3.9)  is  very  efficient  for  viscoplastic  problems  and  for 
viscoelastic  problems  of  the  Maxwell  model  type.  Here  we  use  this  approach 
generally  referred  to  as  the  BIE  method)  for  crack  tip  fields  governed  by  the 
Bodner-Partom  equations.  The  procedure  proceeds  as  follows 

1)  Compute  the  initial  field.  Since  the  inelastic  strain  is 
t  ■  0,  the  initial  field  is  simply  the  linearly  elastic  one.  In  the 
context  this  field  is 

S|t-0  “  %  [KlJl<9)  +  KII?II(0)' 

Hence,  the  solution  is  known  for  t  *  nAt  with  n  =  0,  where  At  is  the 
increment . 

2)  Subdivide  the  near-tip  domain  into  cells,  and  compute 
centroid  of  each  cell  by  using  (5.1).  Then  assuming  that  ei  may  be 
constant  over  each  cell,  u  is  computed  by  using  (3.9)  and  Vu  by 
differentiating  (3.9). 

3)  Compute  x  from  (2.5)  and  (2.6).  Update  x  by  using 

x[(n  +  1) At]  *  x(nAt)  +  At  x(nAt).  (6.2) 

4)  Use  the  updated  value  of  x  to  compute  £•*  as  in  2),  and  repeat  2) 
and  3) . 

When  the  edges  of  the  cells  are  straight,  the  integrals  involved  in 


zero  at 
present 

(6.1) 

time 

at  the 
taken  as 


(3.9)  can  be  evaluated  analytically,  as  discussed  by  Cruse  and  Polch  [11]. 

By  picking  up  the  square-root  singular  terms  in  these  integrals  and  summing 
up  their  coefficients,  ^  can  be  computed.  Alternatively,  the  computation 
of  ^  may  be  carried  out  by  using  (4.27)  with  the  piecewise  constant 
approximation  for 


.7.  Numerical  Results 


Numerical  computations  have  been  carried  out  for  the  non-hardening  case 
of  the  Bodner-Partom  model.  With  reference  to  Eqs.(5.1)  and  (5.4),  this 


case  corresponds  to 

Z  /Z  -  1 
o  1 


(7.1) 


It  is  then  convenient  to  introduce  non-dimensional  quantities  as 
follows 


U  -  u  (yZ//3)/(Kp  ,  T  -  i  /(Z//3) 


(7.2a,b) 


x  -  *  (Z2/3)/(K^)  ,  t  -  t  (uDq)/(Z//3) 


(7.3a,b) 


E  -  —  /D  (7.4) 

dt  o 

Here  the  dot  notation  represents  the  derivative  with  respect  to  t.  In 
terms  of  these  dimensionless  quantities,  Eqs.(5.1),  (2.6)  and  (2.8)  convert 


t  -  [T,/(k,-r,),w-  ?<lr-r>~n] 


T  -  — (7*y)l  +  vu  +  (vu)T  -  2E1 


(7.5) 


(7.6) 


V2U  +  7V-U  -2V-E1  ,  (7 

where  V  and  V2  are  the  del  operator  and  the  Laplacian  with  respect  to  the 
dimensionless  optical  coordinates  x^,  and  T'  is  the  dimensionless  stress 
deviator.  Since  X/v  and  (X+u)/u  can  be  expressed  in  terms  of  the  elastic 
Poisson's  ratio,  v  ,  it  now  follows  that  numerical  values  have  to  be  chosen 
only  for  v  and  n.  Here  we  choose 


(7.7) 


a; 


.8 


2 


For  the  numerical  procedure  based  on  the  boundary  integral  equation 
method  described  in  Section  6,  the  cell  arrangement  shown  in  Fig.  4  has 
been  used.  The  expression  for  the  time  rate  of  change  of  the  stress  intensity 
factor  is  given  by  Eq.(4.27).  Once  the  inelastic  strain  rate  has  been 
computed,  the  term  d*  can  be  determined  and  (4.27)  can  be  evaluated 

numerically.  A  subsequent  numerical  integration  with  respect  to  time 

~  E  - 

produces  K^(t).  Figure  5  shows  K^(t)/K^  versus  the  dimensionless  time  t. 

-  £ 

It  is  noted  that  K^(t)/K^  decreases  with  time.  Also  shown  in  Fig.  5  is 

the  normalized  stress  (x^,0,t)  /T^  (x^,0,0)  ,  evaluated  at  x^=  0.0014.  The 

two  curves  diverge  rapidly  as  the  dimensionless  time  increases. 

It  is  instructive  to  compute  the  actual  physical  distance  to  the  crack 

tip  which  would  correspond  to  x^=»  0.0014.  Suppose  we  consider  a  crack  in  an 

unbounded  medium  under  far-field  uniform  tension  T_„  ■  o  .  Then 

12  o 

g 

Ki  *  ao/(ira).  For  titanium,  Ref.fl]  suggests  the  following  value  for  the  yield 

stress,  aY  ~  250  N/mm2,  while  Z  ~  1200  N/mm2.  Let 

-  a^/2.  Equation  (7.3a)  then  yields  x^~  0.102ax^.  Hence  x^=  0.0014 

-4 

corresponds  to  x^~  1.428x10  a,  which  is  a  very  short  distance  to  the  crack 
tip.  The  two  curves  in  Fig.  5  thus  imply  that  even  very  close  to  the  crack 
tip,  the  singular  term  by  itself  would  give  a  very  poor  approximation  to  the  near- 
tip  stress  an  error  which  increases  with  time  t.  The  relation 

between  physical  time  t  and  dimensionless  time  t  can  be  obtained  from  (7.3b) 
as  t  *.635xl06t,  where  the  following  values  from  Ref.[l]  have  been  used 
Dq  -  lQVsec  and  y  =  . 44x10 SN/mm2 .  As  predicted  by  the  constitutive  model 


used  in  this  paper*  the  relaxation  of  the  stress  intensity  factor  takes 
place  very  rapidly. 

An  improved  asymptotic  representation  of  the  near-tip  stresses  can  be 

obtained  by  Including  higher  order  terms.  The  calculation  of  these  terms 

requires  the  evaluation  of  the  integrals  for  L^,  Eq.(4.28),  and  f_(9  ), 

Eq.(4.32).  The  r  logr  term  in  Eq.(4.23)  contributes  only  a  constant  term 

to  A  constant  contribution  is  also  obtained  from  the  term  f_(0  )r  . 

zz  o  o 

After  the  calculations  have  been  completed  we  find  on  the  crack  line  (6  *  0) 


3.376 


(7.9) 


or 

T22  -  -  3.376t  .  (7.10) 

This  two  term  expansion  gives  excellent  agreement  with  the  numerical  result, 

as  can  be  seen  from  the  numerical  values  listed  in  Table  1. 

_  e  — 

Figure  6  shows  Kjjj(t)/K^j.  versus  the  dimensionless  time  t.  Also 

shown  is  the  ratio  T^x^.O.^/T^Cx^C.O)  evaluated  at  0.0047.  As  for  the  Mode  I 

case,  the  two  curves  diverge  rapidly  as  t  increases. 

It  is  relatively  simple  to  write  out  explicitly  the  asymptotic  results 

for  the  Mode  III  case.  At  the  crack  tip  we  obtain  from  Eq.(7.5) 


^(O.e) 


/  ~a±r¥  \ 
\  cos^e  / 


(7.11) 


Use  of  this  result  in  (4. 36)-(4. 38) ,  and  subsequent  Integration  with 


respect  to  time  yields 


/r 


(7 


S±t*2  \ 

j  siiry  + 

3 

sinj 

0  J 

^  30 

0 

COSJ 

'  — COSTJ - 

3 

COSJ 

The  values  of  at  x^  ■  0.0047,  ■  0  computed  according  to  Eq.(7.12)  and 

by  the  B1EM  have  been  listed  in  Table  1.  Good  agreement  is  observed. 

The  results  presented  in  this  section  show  that  the  square  root  singular 
term  presents  a  very  poor  approximation  to  the  crack  tip  field  even  at  very 
small  distances  from  the  crack  tip.  The  approximation  becomes  worse  as 
time  increases.  The  use  of  additional  terms  in  the  asymptotic  expression 
for  the  particle  velocity  provides,  however,  a  very  satisfactory  agreement 
over  the  whole  range  of  time  between  the  corresponding  asymptotic  represen¬ 
tation  and  the  numerically  computed  stress  component  at  a  crack-line  point 
just  ahead  of  the  crack  tip. 
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Appendix  A;  Finite  Parts  of  Integrals 
The  "finite  part"  of  an  integral,  denoted  by  pf,  is  defined  as 
(see  [14],  pp.  75-104) 

pf  I  f(x)dx  -  lim  [  I  f(x)dx  -  divergent  part  as  R-+®,  e-*o] 

•h  P-MS  Jp 


A  typical  example  is 


r 


riogro  ,  n  -  1 


,  n  >  1 


For  small  r  it  then  follows  that 
o 


f o  o  (r,9) 

pf  J  - n“  dr 


(logro)g  (0,9)  +  o(rQ)  ,  n  -  1 


(-TT  ri’n)oi(O,0)  +  n  >  1 

1— n  o  o 


where  we  have  used  that  d  (r,0)  is  bounded  at  r-0,  see  (4.3).  Consistent  with 

the  definition  (A.l),  it  is  convenient  to  rewrite  (4.18)  in  the  form 

rrfl  r  n+1 

.  JT  08  r  r  ~~T~  *  C°  r  T" 

1>(0,9o)  “  f  d9  (  1  [Pf  /  (T2)  Nn+1(0,0o)o1(r,0)dr  -  pf  J  (^) 

~n  n*l  o  '  o 


-n+l^’V?  (r,0)dr]} 


By  using  (A. 3),  Eq.(A.4)  can  be  rewritten  as 


I>(0»0r>)  ”  f9(®  )*  1°8  r  +  r  B(0  )  +  o(r  ) 

>  o  2oo  o  o  o  o 


where 


1 


has  been  defined  in  (4.22)  and 


B(eo) 


®  g*(r,0) 


CH2(e.Q  )  pf  [  ^ -  dr  +  I  -^■Nn+1(8)8  ).i1(O,0)l  (A.6) 

Jo  n*2 


Equation  (A. 5)  shows  that  the  order  of  I  as  r-K)  is  defined  by  terms  r  log  r 

_>  o  O  C 


and  r  . 
o 


It  is  too  complicated  to  compute  all  the  terms  in  Eq.(A.6)  by  using  the 
expansion  given  by  (4.13).  Another  way  of  computing  the  relevant  terms  follows 
by  considering  Eq.(4.11)  for  the  case  ^(r.O.t)  =  1.  From  (4.11),  (4.14), 
(4.17)  and  (A. 5)  it  then  follows  that 

pf  f  ■  V  I  BlW-vi.,  Ail.,  <•.»_>] 


n— 1  n-1  -n+1  o 
nj*  1 

-  N2(9»9o)  ro  log  ro 


(A.7) 


Next  we  define 


00 

P(0,0  )  »  —  pf  f  rVGdr  +  N„(0,0  )  log  r 
o  r  J  -2  o  6  c 


(A.  8) 


It  can  then  be  verified  that 


4  .  /(<+l)cos0 

1*  /  A  Jft  a  a-*  I 


f,(0  )  -  f  P(e,e  )-a(o,e)de  +  c 

-3  °  Jv  -  o  - 


1 '(<-3)sin0  ' 


(A.  9) 


-COS0 


The  terms  with  the  constants  and  come  from  the  first  integral  of 

•  •  %  * 

(A.6).  The  constant  C0(C,)  vanishes  for  b  having  Mode  I  (II)  symmetry. 
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The  manipulations  that  are  required  to  evaluate  (A. 8)  by  the  use  of  (3.6) 
are  lengthy  and  complicated.  For  the  Mode-I  case  the  result  is  given  by 
Eq. (4.32) . 

Now  that  the  third  term  in  the  expansion  (4.23)  is  known  explicitly, 
namely , 

u2  -  f2(0)rlogr  (A. 10) 

and  it  is  known  that  u^  is  of  the  form 

u3  -  f3(9)r  ,  (A. 11) 

an  alternative  way  of  computing  f3(0)  is  provided  by  returning  to  the 
substitution  method  outlined  by  Eqs. (4.1)-(4.6).  By  collecting  terms  of 
order  r  we  find 

A*[rf3(0)l  -  A*u2  +  V-o*  (A. 12) 

A  direct  calculation  shows  that  rA*u_  is  a  function  of  0  only.  Hence  the 

x-2 

terms  on  the  right  hand  side  of  (A. 12)  are  r  ^  times  a  function  of  0  .  The 
standard  Fredholm  argument  shows  that  for  boundary  conditions  appropriate  to 
Mode  I,  (A. 12)  does  have  a  solution  which  is  unique  to  within  a  term  which 
would  give  rise  to  a  uniform  contribution  to  r^.  The  solution  to  (A.  12) 
can  now  be  determined  by  well  established  methods.  The  result  is  given  by 
Eq. (4.32) . 
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Inelastic  strain  rates  for  uniaxial  tension,  o,  and  various  values 
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